Magnetic field evolution of the quasiparticle interference in a d-wave superconductor 
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Quasiparticle interference in a d-wave superconductor with weak disorder produces distinctive 
peaks in the Fourier-transformed local density of states measured by scanning tunneling spec- 
troscopy. We predict that amplitudes of these peaks can be enhanced or suppressed by applied 
magnetic field according to a very specific pattern governed by the symmetry of the superconduct- 
ing order parameter. This calculated pattern agrees with the recent experimental measurement and 
suggests that the technique could be useful for probing the underlying normal state at high fields. 



There now remains little doubt that hole-doped 
cuprate superconductors below their critical temperature 
T c form a rather conventional BCS superconducting (SC) 
state characterized by a spin-singlet d-wave order param- 
eter. The key mystery in the field is the nature of the 
state that occurs when superconductivity is suppressed 
by undcrdoping, magnetic field, or temperature Re- 
cently observed quantum oscillation phenomena in high 
magnetic fields [2| indicate a metallic state with small 
Fermi pockets which appear incompatible with the stan- 
dard band structure calculations. These results are also 
difficult to reconcile with the angle-resolved photoemis- 
sion (ARPES) data @, i], which instead imply "Fermi 
arcs" , i.e. disconnected segments of a Fermi surface that 
appear above T c close to the nodal points of the d-wave 
order parameter. 

It is possible, in principle, that the normal state 
reached by the applied magnetic field is different from 
the state above T c . However, since ARPES cannot be 
performed in high magnetic fields and quantum oscilla- 
tions are difficult to detect at elevated temperatures, an- 
other technique is needed to settle this puzzle. A good 
candidate is the scanning tunneling probe which can be 
applied at finite temperature as well as in the presence of 
magnetic fields jfjj]. The recently perfected technique of 
Fourier-transform scanning tunneling spectroscopy (FT- 
STS) allows extracting the Fermi surface from the dis- 
persion of the quasiparticle interference peaks 0, 0, 0] 
and yields results in good agreement with ARPES. 

In this Communication we take a first step in this di- 
rection by studying theoretically the effect of relatively 
weak magnetic fields (such that the sample remains in 
the superconducting state) on the quasiparticle interfer- 
ence patterns observed in FT-STS. We find that the field 
causes strong enhancement of a subset of the interference 
peaks that are observed in zero field. The pattern of en- 
hancement, illustrated in Fig. [TJ is closely related to the 
d-wave symmetry of the superconducting order parame- 
ter and agrees with the recent measurements performed 
by Hanaguri et al. 0]. This agreement exemplifies our 
level of understanding the the quasiparticle dynamics in 
cuprates and lays foundation for the future studies in 
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FIG. 1: (Color online) a) Contours of constant quasiparti- 
cle energy ui = O.lf in a d-wave superconductor in the first 
Brillouin zone. We use standard tight binding lattice model 
[!H with t' = -0.3t, A = 0.2f and ^ = -t. b) Positions of 
the quasiparticle interference peaks resulting from the octet 
model downfolded to the first Brillouin zone. Peaks enhanced 
(suppressed) by the magnetic field are marked by solid (open) 
circles. 



much stronger fields. 

In our subsequent discussion we focus on the quantity 
Z{k,ui) measured recently by Hanaguri et al. 0, [lCfl 
defined as the spatial Fourier transform of the ratio 



(1) 



where g(r,u>) is the tunneling conductance dl/dV mea- 
sured at point r of the sample at bias Hlu. A key advan- 
tage of considering the ratio Z(r, to) is that the unknown 
tunneling matrix element connecting g(r,u) to the local 
density of states n(r,u) drops out (provided that it is a 
slowly varying function of w) leaving behind the ratio of 
the local density of states which contains information on 
the intrinsic electronic state of the system. 

There are two very interesting aspects of the above 
measurements Q : (i) The quasiparticle interference pat- 
terns in Z{k,uj) are even clearer and more striking than 
those observed in g(k, uj) in the same sample, and (ii) the 
patterns in Z(k, oS) are sensitive to the applied uniform 
magnetic field in the range to 10T. More specifically, 
with the increasing field intensities of various interfer- 
ence peaks vary in a very specific way. In what follows 
we formulate a theory of this field-induced variation. 
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The local density of states (LDOS) in a material can 
be decomposed into two parts 



n(r, lu) = uq(lu) + Sn(r, lu), 



(2) 



The first part is uniform in space, reflecting the physics 
of a perfectly homogeneous native material and, for a 
d-wave SC is a V-shaped function of co. The second 
part describes inhomogeneity due to disorder. As dis- 
cussed extensively in the literature, the structure of the 
quasiparticle excitation spectra together with the BCS 
coherence factors cause the Fourier transform of 6n(r, lu) 
(which we refer to hereafter as FT-LDOS) to comprise 
a collection of sharp peaks [ll|, [l2|, EH • The location of 
these peaks and their dispersion as a function of lu can be 
understood from a heuristic octet model Q based on a 
set of eight points in the Brillouin zone at the tips of the 
banana-shaped contours of constant energy illustrated in 
Fig.QJa). A peak in FT-LDOS will appear at momentum 
qi if it connects any two of the octet points. 

To gain theoretical insight into the structure of Z(r, lu) 
we now substitute Eq. ([2]) into ^ and expand to leading 
order in 5n 



Z(r,u) ~Z (lu) 
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5n(r,Lu) 5n{r,—uu) 



u (lu) 



rio(-w) 



(3) 



where Z (lu) = n^w) /n${—Lu). Eq. ^ should be an ex- 
cellent approximation as long as |<5ro(r,a>)| <C |n, (w)|, 
a condition well satisfied for the data under considera- 
tion [Io| . We are interested in the spatially varying part 
5Z{r,uj) of the above expression. It is useful to recast it 
in the following way: 

SZ(r,w) = Cx(w)8n e {r,u) + C 2 {uu)8n {r , lu) , (4) 

where Sn ^ represents the part of Sn even (odd) in lu 
and Ci(2) are r-independent functions of frequency. 

The above even/odd decomposition facilitates the fol- 
lowing key observation. As pointed out by Chen et al. 
14 1, for a strictly particle-hole symmetric system, 5n a 
originates exclusively from the scattering in the particle- 
hole channel (i.e. ordinary potential scattering) while 5n c 
comes from scattering in the particle-particle channel (i.e. 
modulation in the SC order parameter). Cuprates are of 
course not strictly p-h symmetric but nevertheless they 
are sufficiently close to the p-h symmetric situation that 
the above classification holds to a very good approxima- 
tion In the absence of the magnetic field disorder 
in the sample gives rise to ordinary potential scattering 
as well as off-diagonal scattering caused by local suppres- 
sion of the SC gap amplitude by pair-breaking impurities. 
Thus, Sn has in general both even and odd contributions, 
even in the strictly p-h symmetric case. 

When the magnetic field is applied to the sample in ex- 
cess of the lower critical field H c \ the sample enters the 
mixed state and Abrikosov vortices appear. If there is 



significant randomness in the Abrikosov lattice then vor- 
tices cause additional quasiparticle scattering due to (i) 
the suppression of the order parameter in the vortex core 
and (ii) the superflow associated with the screening cur- 
rents outside the cores. These effects both contribute to 
scattering in the particle-particle channel and thus pre- 
dominantly affect Sn c . Ordinary potential scattering, by 
contrast, should be largely unaffected by the magnetic 
field. We thus expect magnetic field to enhance Sn e but 
leave Sn Q largely unchanged. 

In the following we shall explicitly calculate 5n(r,u>) 
caused by the order parameter suppression in the vortex 
core, which we believe is the dominant effect of the mag- 
netic field. We verify that it is indeed predominantly 
even in frequency in the vicinity of the p-h symmetric 
point and analyze in detail the spatial structure of the 
resulting interference pattern reflected in Z(k,u>). 

The local density of states in a superconductor is given 
by n(r, lu) = — -lmGu(r, r; lu + id) where G(r, r'; lu) is 
the full electron Green's function, a 2 x 2 matrix in the 
Nambu-Gorkov space. For weak impurity scattering the 
FT-LDOS can be calculated using the Born approxima- 
tion [H 

Sn(q) = --V a (q)ImJ2{G (k)<j a G°(k - g)] Uj (5) 

k 

where o~ a are Pauli matrices in the Nambu space, V a (q) 
is the Fourier transform of the random impurity potential 
in the charge (a = 3) and spin (a = 0) channel. The k 
summation extends over the first Brillouin zone and the 
frequency arguments have been suppressed for brevity. 
G° denotes the unperturbed Green's function 



G (k,iw) = 
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El 



ilu - e k 



(6) 



with £fe the band energy measured from the Fermi sur- 
face, Afe the gap function and Ek = \AI + A|. 

The LDOS modulations Sn(k, lu) due to im pur ity scat- 
tering have been extensively studied [ll, 12, 13, 15, li| 
based on Eq. ([5]) as well as more accurate t-matrix cal- 
culations. Comparison to a series of atomic resolution 
FT-STS data [aS IH shows good qualitative agreement 
in terms of peak positions and dispersions. 

When considering scattering off of spatial modulations 
of the SC order parameter such as those occurring near 
the vortex core, one might expect that a formula just like 
Eq. ([5]) but with off-diagonal Pauli matrices (a = 1,2) 
should be applicable. This would indeed be the case for a 
simple s-wave superconductor. In the case of a <i-wave or- 
der parameter the situation is slightly more complicated 
[Tit - This is related to the fact that the d-wave order 
parameter is most naturally described as living on the 
bonds of the underlying square lattice. Correspondingly, 
a point-like perturbation will be a gap modulation SAi 
that affects four bonds emanating from a single site r^. 
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A general gap modulation can be thought of as a sum of 
these point-like modulations. 

To formulate this we consider a perturbation described 
by the Hamiltonian 



1 



SH = - Y^SA i xs[c^(r i )c l (r i +S)-c l (r i )c^(r i +S)+h.c.) 
1 i.S 

(7) 

where c a (r) represent the electron annihilation opera- 
tors, 8 = ±x,iy, and xs = 1 f° r z-bonds and —1 
for y-bonds. For simplicity we consider <5A^ real. If 
we define the usual Nambu spinor operator ip( r i) = 
[c T (ri),c|(r 4 )] T we may write SH = J2kM> ^kYkk'^k' 
with 



V k k' = eri<5A fc _ fc /(xfc +Xk>) 



(8) 



and Xk = cosk x — cosfc. y the Fourier transform of xs- 
Had we allowed SAi to have imaginary part there would 
be an additional component of V kk ' proportional to a- 2 . 

Within the Born approximation we thus have the fol- 
lowing vortex-induced LDOS modulation 



Sn(q) = --Tm V[G°(fc)y fc , fc „ q G°(fe - q)] n . 

7T ' * 



(9) 



If we now identify 5A q with V\(q) we see that indeed 
this result has the form of Eq. ([5]) , except for the factor 
(\k + Xk-q) implied by Eq. ((8]) that reflects the d-wave 
symmetry of the order parameter. Using Eqs. ((6|) and ([8|) 
we may further rewrite 5n(q,ui) = — ^6A q ImA(q, uj + iS) 
with 



A(g,iw) = 5^(x+ + X-) — 

k 



+ e + )A_ + (zlj + 6_)A_| 
(uj 2 + E 2 + ){uj 2 + W) 



(10) 

and e± = e k ± q/2 etc. 

In the particle-hole symmetric limit, the terms with 
e± vanish upon the momentum summation. To see this 
note that in this limit the band energy has the property 
(k+Q = —tk for Q = (it, it). Since Xk and A*, also share 
this property (irrespective of p-h symmetry) it follows 
that shifting the summation variable by Q reverses the 
sign of e± , which therefore vanish in the sum. In the p-h 
symmetric case we are thus left with 



A(q,iuj) 



1 

A~ 



k 



iu(A + + A_) 2 
[uj 2 + E 2 + )(uj 2 +Eiy 



(11) 



where we expressed the gap function as A^ = A \k- As 
will be evident shortly, the remaining expression is even 
in uj when analytically continued to real frequencies. 

We now wish to examine the effect of this contribu- 
tion on the quasiparticle interference peaks that are pre- 
dicted by the octet model. To this end it is useful to 
perform analytical continuation and explicitly evaluate 



— ilmA(g, uj + id) at wavectors 
are the octet vectors. This yields 



Qi " Qj ) where Qi 



1 



k 



8{\u\-E k )V 



(Afc_ g + Afc) 2 

(El- 



El) 



(12) 



where V denotes the principal part. It is clear that for 
q = qtj the largest contribution to the sum comes from 
the vicinity of k = Qi, Qj. Near these points the denom- 
inator approaches zero but numerator is a slowly varying 
function of k. We may thus approximate the latter by 



its value at k = 
We thus obtain 

- -ImA(qry) 

7T 



Qi , Qj and take it outside of the sum. 



(At + Aj) 2 — 



V 



5(\uj\-E k ) 



(13) 

where Aj denotes A q evaluated at the octet vector Qi. 

The above Eq. ()13|) has some remarkable implications 
and represents our main result. Most importantly it im- 
plies that in the p-h symmetric case the effect of the 
applied magnetic field on the octet vectors can be sum- 
marized as 



5n c {q v 



{Ai + AjfKijiuj). 



(14) 



Here /Cy (uj) denotes the sum in Eq. (| 13|) and can be 
shown to represent a positive quantity whose precise 
value for a given frequency and vector qfy depends on 
the details of the underlying band structure and is thus 
non-universal. The factor (Ai + Aj) 2 is, by contrast, uni- 
versal and depends only on the symmetry of the SC order 
parameter. Since the octet points lie on the Fermi sur- 
face of the underlying normal metal it is easy to see that 
A; | = \uj\ for all Vs. The sign of Aj, however depends 
on the position of Qi in the Brillouin zone as illustrated 
in Fig.QJa). It then follows that 



8n e (qij,uj) 



Auj 2 JC tj (uj), if sgn Ai = sgn Aj , , 
0, if sgn Aj ^ sgn Aj 



Thus, remarkably, we find that only those interference 
peaks will be enhanced by applied magnetic field whose 
wavectors qij connect octet points in the regions of the 
Brillouin zone with the same sign of the gap function A& , 
denoted by +/+ and — /— in Fig. Q^a). These are q\, q4, 
and Q5. The remaining +/— peaks will be to leading 
order unaffected. The resulting pattern is illustrated in 
Fig. \V[b) . We remark that these are precisely the peaks 
observed to be enhanced in the experiments by Hanaguri 
et al. [§]. 

If the system breaks the p-h symmetry, as is generally 
the case in cuprates at finite doping concentration, then 
the above conclusion cannot be formulated as a precise 
symmetry statement. However, as long as the p-h sym- 
metry breaking remains relatively weak, as is the the case 
in cuprates close to half filling, our result Eq. (fT"5)) will 
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FIG. 2: (Color online) Enhancement of FT-LDOS peaks by 
magnetic field at octet vectors q; modeled by Eq. (|10|l with 
the band structure as in Fig. [T] In the calculation the field H 
is represented by the vortex core scattering rate Vi and Ho is 
chosen such that H /H0—V1/V3. The data are normalized to 
the H = value which we model by Vi = and V3 > 0, i.e. 
charged impurities only. 



hold to a good approximation: the +/+ and — /— peaks 
will be significantly enhanced while +/— will be affected 
only slightly. This can be seen by analyzing Eq. (JTUJ). 
Even when p-h symmetry is weakly broken we can still 
write it as (fTTj) plus a correction that will be odd in fre- 
quency. This correction will be (i) small compared to the 
leading term by factors of fi/t and t'/t, where fj. is the 
chemical potential and t, t' are nearest and next-nearest 
neighbor hopping amplitudes, and (ii) will contain a fac- 
tor (A + + A_) which will make it very small for the +/— 
peaks, as before. 

In order to ascertain the validity of our above conclu- 
sions we have evaluated the sums indicated in Eq. (fT0|) 
numerically for band structures with realistic parame- 
ters. This is illustrated in Fig. O We have also analyzed 
these expressions within the nodal approximation along 
the lines of Refs. [l3|, [l5|. These considerations confirm 
that Eq. (fTSj) is an excellent proxy for Eq. (jTUJ) when the 



p-h symmetry is present as well as when it is weakly vio- 
lated, as in the superconducting state of cupratcs. Specif- 
ically, we find that magnetic field enhances the peaks at 
tjfi, q.4, and q$, while the remaining peaks are essentially 
unaffected. This pattern of enhancement remains sur- 
prisingly robust even when the p-h symmetry is signifi- 
cantly violated as in Fig. O 

Hanaguri et al. Q report that the +/+ and — /— are 
enhanced by the applied magnetic field and the +/— 
peaks arc in fact reduced in amplitude. This can be 
reconciled with our theoretical prediction when we re- 
call that our model explicitly treats only one aspect of 
the field, namely the suppression of the order parameter 
in the vortex cores. Magnetic field also generates su- 
perflow which is known to produce a Doppler shift (l7l | 
and a more subtle Berry phase [HI effect on the quasi- 
particle wavefunctions. Since these are both long-range, 
non-local effects, their impact on the quasiparticle in- 



terference patterns is significantly more difficult to com- 
pute. It appears to us likely, however, that the additional 
phases acquired by the quasiparticles as they propagate 
on the background of the random vortex array will tend 
to scramble the interference patterns and therefore sup- 
press the peaks. We thus hypothesize that a combination 
of this suppression of all peaks and the enhancement of 
+/+ and — /— peaks due to the vortex core scattering 
will lead to the pattern observed in experiment Q . 

Our results here underscore once again the importance 
of the quasiparticle coherence factors [HI for the tunnel- 
ing interference spectroscopy. Indeed, the pattern of the 
peak enhancement by magnetic field found here is deter- 
mined solely by the coherence factors. Their presence, 
manifested in the peak-like FT-STS patterns, indicates 
pairing even in magnetic field. Our results, in conjunc- 
tion with experimental data Q, also illustrate the re- 
markable sensitivity of the FT-STS technique to rela- 
tively modest magnetic fields up to 10T. This suggests 
good prospects of FT-STS for unraveling the mystery of 
the normal state reached at higher fields or temperatures. 

The authors are indebted to J.C. Davis and T. 
Hanaguri for discussing their data prior to publication 
and to P. Coleman, E.-A. Kim and L. Taillcfcr for illumi- 
nating discussions. This work was supported by NSERC, 
CIFAR and the Killam Foundation. 
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